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We present an analysis of the magnetic susceptibility and specific heat data of Kapellasite 
(Cu3Zn(OH)eCl2) through high-temperature expansions (HTE) of the J\- J2- Jd -Heisenberg model 
on the Kagome lattice (up to order 9). Experimental data are reproduced with exchange energies 
Ji = -12.4(2) K and J 2 = -5.0(8) K and J d = +14.8(8) K, where J d is the third neighbor exchange 
energy across the hexagon. The agreement between theory and experiment is excellent for the spin 
susceptibility down to 20K. Analysis of the spin specific heat data shows a deficit of 10% of the 
entropy around 10K. This improved modeling confirms the results of ref 3 regarding the location of 
Kapellasite in the cuboc2 phase of the Heisenberg model. The quality and limits of this modeling 
are discussed. 

PACS numbers: 02.60.Ed 71.70.Gm 75.10.Kt 75.30.Et 



I. INTRODUCTION 

Kapellasite first synthesized in London 1 2 is a poly- 
morphous of Herbertsmithite of chemical formula 
Cu3Zn(OH)6Cl2. Contrary to Herbertsmithite, the sus- 
ceptibility of this new metastable compound points to a 
ferromagnetic Curie Weiss field of about 10 K. Kapell- 
asite was not found to develop dominant ferromagnetic 
correlations down to the lowest temperature and this in- 
formation is the landmark of competing interactions. As 
this compound fails to develop any order down to 20 
mKP a quantitative discussion of the possible spin liq- 
uid properties needs a modeling of the spin hamiltonian. 
Such a model could be built by successive approxima- 
tions from a fit of the thermodynamic quantities to high 
temperature series expansions for a relevant set of hamil- 
tonians. The actual chemical formula of the synthesized 
compound, determined with neutron powder diffractiorP, 
is (Cu .73Zno.27)3(Zno.88Cu .i2)(OH) 6 Cl 2 , with 27% Zn 
on the Cu-sites of the kagome lattice and 12% Cu on the 
hexagonal Zn site. The Cu/Zn mixing leads to some dis- 
order within the kagome planes but cannot induce any 
coupling between the planes, which occurs only via very 
weak O-H-Cl hydrogen bonds^. Kapellasite is therefore 
remarkably two-dimensional. 

A first theoretical description of Kapellasite, which is 
deep in the Mott phase, is the Heisenberg Hamiltonian 
on the perfect kagome lattice: 

H=J2 J«Si-Sj, (1) 

where the exchange integrals J a are defined in Fig. [T] 
Due to the geometry of the exchange paths J3 and Jd are 
different and Jd is expected to be larger than J3 by an 
order of magnitude 5 . We will thus limit this first analysis 
to a pure J\-J2~Jd -model, neglecting both the effects of 
disorder and of an eventual Dzyaloshinskii-Moriya inter- 
action. These issues will be discussed in the conclusion. 




FIG. 1. (Color online) Kagome plane of Kapellasite with 
Cu 2+ S— 1/2 spins (blue), non-magnetic Zn 2+ ion (green), 
and exchange interactions (red). 



Because of the competition between these various ex- 
changes, the shape of the susceptibility and heat capac- 
ity versus T curves deviates noticeably from those of the 
pure Heisenberg model. In fact the 17th order HT series 
expansions of the nearest neighbor Heisenberg hamilto- 
nian 6 deviates significantly from the kapellasite data for 
temperatures lower than 70 K. We thus compute the HT 
series of magnetic susceptibility X and specific heat C v 
with the J\-J2~Jd parameters up to order 9 (see appendix 
A). 

The paper is organized as follows. In section II, the 
magnetic susceptibility X is fitted to experimental data 
providing strong constraints on the coupling constants. 
Section III is devoted to the specific heat analysis allow- 
ing us to further refine these constraints. Section II and 
III are organized similarly. We first detail the technical- 
ities of the fits and finish with the physical conclusions 
to be kept in mind for the appreciation of the physical 
properties of the model. In the conclusion we discuss 
the consequences of neglecting at this stage the chemical 
disorder in the plane and Dzyaloshinskii-Moriya interac- 
tions. 
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FIG. 2. (a) and (b): Fit quality, as denned by 
Eq|6] with e = 10 -3 , versus V2 (a) or z/3 (b) for 
T min = 28.5 K (a) and T min = 19.5 K (b). Contour 
levels are at integer values. The best fit is for the high- 
est value of Qx (thus in red). The best points are 
{Ji, u 2 , A, B}={-19.34K, -0.483, 1.012, -6.5 x 10" 5 K' 1 } 
(a), {Ji, z/ 3 , A, B}={-U.74K, -0.850, 1.023, -9.1 x 
lO -5 ^ -1 } (b). (c) and (d): comparison with experiment for 
these best points of (a) and (b) respectively. All PPAs at 
order n = 10 (c) n = 11 (d) are shown and d indicates the 
degree of the denominator for each PPA (see Appendix [b] for 
the PPA's definition). The thick vertical line indicates T m i n . 



II. DESCRIBING X(T) 

The dc susceptibility was measured in a commer- 
cial Quantum Design MPMS-5S SQUID magnetometer. 
Whereas a divergent behavior of x is commonly ob- 
served in most frustrated magnets and interpreted as 
a signature of weakly coupled magnetic defects^, here, 
the NMR local probe data indicate that the measured 
macroscopic susceptibility is intrinsic. The experimental 
data are given as a list of points {T&, A^ xp }. Fitting the 
data to a Curie- Weiss law X{T) ~ C/{T - T ), leads to 
C = 0.429(2) cm 3 K/mol and T = 9.5 ± 1 K, where T is 
interpreted as the Curie- Weiss temperature. In the range 
of temperature of interest XT /C is of order of unity and 
is thus a good quantity to fit. 

We define the HT-series expansion of the magnetic sus- 
ceptibility X UT : 



X (T)T 
C 



(2) 



where f3 = 1/T and Pi is an homogeneous polynomial 
of order i and n is the highest order at which the series 



is known. Pi(Ji, J2, Jd) = and for the Kagome lattice 
= — Ji — J2 — Jd/2. Assuming J\ ^ 0, we use as 
parameters = J 2 1 Ji and V3 = Jd/ J\ instead of J 2 and 
Jd, and define the polynomials ^(^2,^3) (see Appendix 
|Al) as: 



<Y H1 (T)T 
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(3) 



In order to account for the uncertainties in the num- 
ber of spins and the temperature independent Van Vleck 
and diamagnetic susceptibilities, we introduce two pa- 
rameters A (close to 1) and B and define a least mean 
square error as: 



k=i 

Rx = v Zx/N T 



i=0 



BT h 



c 



(4) 
(5) 



where Nt is the number of experimental measurements 
with T > 

Jmin? T\ = T m i n is the lowest tempera- 
ture included in the fit and R x is a fit-error measure. 
One could minimize Rx with respect to the parameters 
{A, 5, Ji, V2-, ^3} taking advantage of the two linear pa- 
rameters A and B. But using Pade approximants, allows 
to extend the fits to lower temperatures. 

In the present method the best set of parameters is 
searched amongst that having the largest number Pade 
approximants providing a "good fit", i.e. the lowest Rx 
error for the largest interval of temperatures. We thus 
define a measure of the fit quality as: 



Q 



}x = ^ Qx,PPA 
{PPA} 
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^,PPA 
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(6) 
(7) 



where the sum runs over the Physical Pade Approximants 
(PPA), see appendix |bJ e = 2 x 10 -3 is of the order 
of the experimental uncertainties while the parameter r 
controls how much the values outside of the uncertainties 
e do contribute to Qx- A large value of r enforces the 
weight of "good" PPAs while it discards more rapidly 
the worse ones (here, r = 8). Rx,ppa is denned by re- 
placing the HT-polynomial by a PPA in Rx (see Eqj4|. 
If a PPA goes through all the experimental points, then 
Q^,ppa is of the order of unity. On the other hand, if 
the PPA is "bad", Rx,d* /e is much larger than 1 and the 
contribution to Q^ppa will be small. Thus Qx counts 
roughly the number of "good" PPAs. For T m i n > 2b K 
we find that Qx is very close to n + 1, i.e. each pade 
approximant is a PPA and i^^,ppA < e. A large domain 
of parameters provides similar quality for the fits. This 
domain is restricted by using smaller values of T m - in . 

Unfortunately, a minimization algorithm is not appli- 
cable because the function Qx is not continuous (see ap- 
pendix p| and has many local minima. On the other 
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FIG. 3. (a): Fit quality, as denned by Eq[6]with e = 10~ 3 , 
versus is 3 and zy 2 for T m i n = 19.5 K. The different narrow 
domains correspond to the values of J\ given in parenthe- 
sis. The dashed line separates the ferromagnetic phase from 
the cuboc2 phase of the corresponding classical Heisenberg 
model. Note that the corresponding — — Jl(l + V2 + ^3/2) 
varies from — 8AK to — 10. 4 X from the top- left solution to 
the down-right one. (b) and (c): the twelve sublattice-cuboc2 
classical structure: the spins of a triangle in (b) start from (or 
arrive to) a summit of the tetrahedron in (c). Spins around 
an hexagon are on a facet of the tetrahedron and spins across 
an hexagon ( Jd interaction) are on the same side of the tetra- 
hedron, thus anti-parallel. 



hand, as the number of parameters is reduced, the qual- 
ity function Qx can be evaluated on grids and after some 
trials the main minimum is eventually found. 

At high temperature, almost all PPAs coincide with 
the HT-polynomial. As the temperature decreases, the 
PPAs start to deviate from each other or deviate from 
the experimental data and thus the quality of the fit de- 
creases. As the function Qx(T m [ n ) decreases sharply at 
some T c , in the following T m [ n is chosen just above T c . 

At fixed the polynomials ^(^2,^3) are 
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FIG. 4. (a): Fit quality, as defined by Eq(6]with e = 10" 3 
with T m in = 19.5 K. The linear relations are defined by Eqs|8} 
[9] with the parameters reported in Table [Tl] (b): Comparison 
with experiments where all PPAs are drawn. 
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FIG. 5. Same as Figgwith T min = 13.5K. 



evaluated to numbers and x UT T/C becomes a polyno- 
mial of f3 = 1/T only. Then all PPAs are computed (see 
appendix |B| . The parameters A and B must be kept 
fixed for all PPAs and their calculations are explained in 
Appendix [C] Thus for a set of given coupling constants 
a value of Qx is computed and we can look for the best 
coupling constant set. 

The pure kagome model {Ji,0,0} where the series is 
known at order 17 is compatible with the experimen- 
tal data with ferromagnetic J\ ~ — 12 if, A — 1.037, 
B = -1.2 x 10-^K- 1 but only for T > 70K. Then 
we check the models { Ji, ^2, 0} and { Ji, 0, ^3} where the 
polynomial has be computed up to order 10 and 11 re- 
spectively (see Appendix |A| ). Fig. [2]-(a)-(b) shows Qx 
while Fig. [21(c)- (d) shows all PPAs at the best points 
of (a)-(b). The fits are of similar quality but for model 
{Ji,^2,0} T min is significantly higher. This is not be- 
cause the series is known at higher order but because the 
shape of XT around the maximum is different. Note that 
for these two models J\ is ferromagnetic while V2 (or ^3) 
is negative, thus J 2 (or Jd) are anti- ferromagnetic. In 
each case the precision on Vi is an order of magnitude 
better than for J\. 

The next model involves {^1,^2,^3}- Fig. 3] shows Qx 
for various values of J\ as a function of and v<$ . We first 
see that for each J\ only a narrow region leads to good 
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-9.1 x 10~ 5 i^ _1 


-7.6 x 10 5 if -1 


Ji 


—13.76(15) K 


—17.89(4) K 


V2 


+0.138(20) 


-0.300(2) 


v-z 


— 1.042(25) 


— 0.oo4(3J 


J 2 


—1.9(3) K 


+5.37(4) K 


Jd 


+14.4(5) K 


+6.51(7) K 





8.5(5) K 


-9.3(1) K 


Jc v 


10.5(3) K 


11.78(4) K 



TABLE I. Results of the fits of Fig(2]5] Errors on Li are 
determined from the localization of the maximum to within 
0.01. and Jc v are the coefficients of the leading term of the 
magnetic susceptibility and specific heat at high temperature. 
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FIG. 6. Low temperature behavior of the Kapellasite specific 
heat. Points are experimental data and the line is a linear fit. 



fits. However, it is clear that two combinations of the 
parameters are strongly constrained. For example, we see 
that in each map, the maxima are near a line ^2 + ^3 = 
/(Ji), with /(Ji) « -a Ji - 1.5, with a = 0.05 K' 1 . In 
order to capture the best domain, we change variables 
using an orthogonal transformation: 




M x 




(8) 



where v\ is dimensionless. In this equation M = 7Z x T 
is a 3 x 3 matrix, where T defines a new orthogonal basis 
with L3 = v\ + V2 + and two orthogonal vectors to L3, 
while 7Z is a rotation of these two vectors around L3: 
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(9) 



where <\> is chosen to obtain a domain {^1,^2} as com- 
pact as possible. For fixed L3, Qx is computed for a 
range of values for L\ and L2. Figure [4]- (a) shows Qx 
at the best L3 while [4|(b) shows the comparison with 
the experimental susceptibility at the best parameters 
for T m i n = 19.5 K. Values of the parameters and results 
for the best solutions are given in Table [TTJ If we extend 
the fits to lower temperatures, a second point emerges, 
not as good around the maximum of XT. FigjE] shows 
this second possible solution. With this procedure we 
find two relations between the exchange energies which 
are known with good precision (Li and L3) while a large 
uncertainty remains on the last one (£2)- 

We finish this section with a few remarks: i) the un- 
certainty in the number of spins, measured here by A 
in Eq|4j is about 2% a value in agreement with experi- 
mental estimates. This shows that this model is realis- 
tic, ii) The Van Vleck and diamagnetic susceptibilities, 
B in Eqj4j less than lO -4 ^ -1 , is also of the order of 
experimental uncertainties, in) both and Tq have a 



positive sign usually pointing to a ferromagnetic phase, 
but all the solutions found in this section fall in the anti- 
ferromagnetic Cuboc2 domain of the classical phase di- 
agram of the J\-J2~Jd model; iv) the values of found 
with the high order fit and Tq coming from the simple 
Curie- Weiss law differ slightly, due to neglecting the Van 
Vleck and diamagnetic susceptibilities (B in EqQ in the 
Curie- Weiss law; v) the fit of the spin susceptibility alone 
does not determine uniquely the three parameters of the 
model but constrains them strongly (see Figj3|. 



III. DESCRIBING CV(T) 

Experimental data are given as a list of points 
{Tfc,Cy X £}. At high temperature, the leading term of 
the spin specific heat decreases as J^/T 2 , where J^ v is 
a positive quadratic form of the coupling constant { J a }, 
here J% v = 3/8(J x 2 + J 2 2 + Jj/2). On the other hand, 
at low temperature, the phonon-specific heat starts as 
(T/Td) 3 , where To is a Debye temperature. When Jq v 
is much smaller than Td, as for helium-3^, each term 
can be evaluated separately. A quick analysis of the 
Kapellasite phonon specific heat reveals that Tjj ~ 160 K 
whereas, according to the results found in the previous 
section, Jc v ~ 10 K. Thus, between 10 and 100 if, both 
contributions are of the same order of magnitude. In this 
section, we analyze the spin specific heat at sufficiently 
low temperature such that the phonons can be neglected. 

In the following cy denotes the specific heat per spin. 
The experimental specific heat per spin is: 



N s k R 



(10) 



where N s is the number of spin per mole. The leading 
low-temperature behavior of the specific heat (essentially 
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a spin contribution) is assumed to be: 



c v (T) = c T a + O(T 



(ii) 



Fig-|6] shows polynomial fits with experimental data of 
the Kapellasite, using a — 2 as expected for a two- 
dimensional anti-ferromagnet and compatible with a 
cuboc2 structure found in the previous section. We find 
the extrapolated value at T = of Cy P (T)/T 2 = Co ~ 
0.075(3) K~ 2 . Assuming a smooth behavior at low tem- 
perature no extra-entropy is expected to be found in this 
regime. We now fit the part of the experimental data 
where the phonons are expected to be negligible, that is 
for T < AK which includes the peak of c v (T)/T. 

Similarly to the magnetic susceptibility (Eqj3|, the 
HT-series expansion of the spin specific heat reads: 
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- 6 




J^-13.7 v 2 =0.14 v 3 = 
J 1 =-17.9v 2 =-0.3 v 3 =-0.364 




\T) = 



2=2 



2T 



(12) 



where qi are polynomials of order i (see Appendix |A| . 

For two dimensional spin hamiltonians, the function 
s(e) is more suitable than cy(T), where s and e are 
the entropy and energy per spin 6 . Indeed, using sum 
rules, s(e) is constrained to start at the ground-state en- 
ergy eo with an entropy s = and end at e = and 
s = In 2, at infinite temperature. Moreover this is a 
monotonic increasing function, (3 = 1/T = s'(e), with 
negative curvature, cy = —s'(e) 2 /s"(e). Imposing the 
low temperature power law for cy(T)^^ ex (T/To) a im- 
plies s(e)e^e0 cx (e — eO) 1 /^, with \i = 1 + 1/a, where 
here a = 2. Thus the function G(e) = s(e) M /(e — eo) is 
analytic in the interval [eo,0]. From the HT-series ex- 
pansion of c S y m (T) we deduce the HT-series expansion of 
G(e) around e = as explained in Appendix [P) 

The HT-series expansion of G(e) is converted into pade 
approximant, noted Gd(e) (d being the degree of the de- 
nominator), and only physical ones (PPA), noted Gd* (e), 
are kept, i.e. those with no zero in the numerator and 
denominator in the interval [eo, 0]. It is sufficient to look 
at Gd*(eo) to compare different PPAs^ indeed all pades 
have the same series around e = and if they have the 
same value at eo it is likely that their variations will be 
very similar. For anti-ferromagnetic ground states, eo is 
usually unknown. Thus this parameter is determined by 
maximizing: 



bOo) 



EE- 

d>2>d* 1 



1 



G d *(e )-C7 d *(e ) 



(13) 



where eo = 10 -3 is of the order of the dispersion G^*(eo) 
in the group of close values and r = 8. Note that the 
found value of eo depends only on the HT-series and is 
thus independent of the fit procedure. Unfortunately, 
this function may be discontinuous because the number 
of PPAs may change eventually. Then the maximum of 
Q e (eo) is found after a systematic search on a grid (see 
Fig0. 



FIG. 7. Variation of Q e versus the ground state energy eo 
for the two sets of parameters determined in the previous sec- 
tion. Despite the complicated variations, one clear maximum 
emerges in each case. The number of points in the interval 
[—9.1, —7.7] is 1000, while the figure is unchanged with ten 
times more points. 



For a fixed set of parameters, once the PPAs of G(e) 
and eo are determined, the spin specific heat are calcu- 
lated for each PPA, noted Cyp PA , as described in Ap- 
pendix [D] (EqsjD8][D10|. To compare with experimental 
data we keep only the data at T < AK such that the 
phonon contributions can be neglected. The exchange 
parameters are determined by maximizing the quality 



factor Qc v as follow: 



i 



PPA 



\Rc v ,ppa/£\ 



Zc v ,PPA = D 

ax l- 



Rc v ,ppa = \l Zc v ,ppa/N t 



(14) 

2 

(15) 
(16) 



where Nt is the number of experimental points with Tk 
less than T max and we used here e = 2 x 10 -3 and the 
parameter D, accounting for both mass uncertainty and 
possible missing entropy, is evaluated as explained in Ap- 
pendix [C] 

For the {Ji,^2,0} and {Ji,0, ^3} models, no good fit 
can be found for the specific heat around the domains 
found for the magnetic susceptibility (see Figj2|. 

Computing Qc v when the eo is unknown is much more 
demanding and less stable and the figures Qc v (^^2^3) 
present lots of discontinuities. Thus proceeding as in the 
previous section is hopeless. Therefore, we focus of the 
constrained domain of parameters obtained from the X 
analysis (Figj3|. At fixed J\ and V2 we look for v<$ that 
maximizes Qx + Qc v - Figj8] shows that the global best 
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FIG. 8. (a): Optimized values of Qx + Qc v versus ^2 + ^3 for 
different values of J\ . (b) : z/2 and z/3 versus J\ for the optimal 
Qx + Qc v • The green shaded domain indicates the overall 
optimal values: J\ — —12.4(2), z/3 = —1.19(3), ^2 + ^3 = 
-0.947(7). 
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FIG. 9. Comparison with experiments for J\ — —Y1AK 
v<i — 0.24 and z/3 = —1.1883. (a): Magnetic susceptibility 
with A = 1.028 and B = -1.05 x 10 
with D = 0.87 and e = -15.53684 K where PPA 
slightly different from the others and PPAs-d = and d — 9 
are not shown as being completely off. The vertical line stands 
for T max (see Eqfl5). 



4 . (b): Specific heat 
1 is 



fit is obtained for : 

Ji = —12.4(1) K 
v 3 = -1.20(4) 
v 2 + ^3 = -0.95(1) 
Jd = 14.8(8) K 
J2 + Jd = 11-8(4) K (17) 
J 2 = -3.1(1.2) K 
= 8.1(5) K 
Jc v = 10.1(4) K 

Figj9]shows the best overall comparison with experiments 
at this point. 

IV. CONCLUSION 

We have fitted the spin contribution of the magnetic 
susceptibility and specific heat experimental data with a 
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FIG. 10. Optimal parameters of the fits in the J1-J2- 
Jd classical phase diagram (Ji = —1) of refP. All fits are 
in the cuboc2 phase. The blue dots are the estimations of 
the best two-parameter models. Grey lines represent the nar- 
row domain, at fixed Ji, where Qx > 5 (Eqj6|. The green 
squares represent the best two values obtained from the fit 
of the magnetic susceptibility. The red large dot is the best 
compromise for the magnetic susceptibility and the specific 
heat. Numbers in parenthesis stand for the corresponding J\ 
values. 



spin-1/2 Ji- J2- Jd Heisenberg model on the kagome lat- 
tice (see Figjl]). In contradiction to ab initio calcula- 
tions of Janson 5 the nearest neighbor coupling is fer- 
romagnetic. This is at variance with Herbert smit hit e 
where the nearest neigbbor interaction is strong and anti- 
ferromagnetic. This can be traced back to the Cu-/i3 QH- 
Cu bonding angle being ~ 13° smaller in kapellasit d 2 * 10 *. 
The related compound Haydeeite Cu3Mg(OH)eCl2, has 
also a ferromagnetic first neighbor interaction but is in 
the ferromagnetic domainP^ This is not the case of 
kapellasite where the J2 and Jd exchange couplings com- 
pete to form a non magnetic compound. 

The spin susceptibility is relatively easy to reproduce 
and imposes strong correlations of Ji, J 2 and Jd along a 
narrow line (see FigfTo]). All solutions stay in the anti- 
ferromagnetic cuboc2 domain of the classical phase dia- 
gram (we do no expect the ferro-cuboc2 transition line 
to move significantly for spin- 1/2 spinsp. The main dis- 
tinctive features of the specific heat data are the low- 
T downturn in C v /T at about 2K, characteristic of the 
competitive exchange couplings and a clear T 2 depen- 
dency excluding a ferromagnetic ground state. The peak 
strongly constrains the parameters. The best domain 
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for both the magnetic susceptibility and the specific heat 
corresponds to J\ between -12 and -14 K, J2 small and 
ferromagnetic and a large anti-ferromagnetic (see red 
dot in Fig 10). Thus competitive exchange energies give 



here a ferromagnetic behavior of the magnetic suscepti- 
bility at high temperature but an anti-ferromagnetic one 
at low temperature for specific heat. 

However the agreement between experiment and the- 
ory is not as good as a quick glance on Fig. [9] would sug- 
gest. There is about 10% of missing entropy (D = 0.87 
in Eqjl5] whereas the mass uncertainties is of only about 
a few percent) in our description. Considering improba- 
ble to find this missing entropy at very low temperature, 
we have to find it at intermediate temperature between 
2 and 20 K where we have not succeeded to fit the full 
specific heat variations with this spin model and phonon 
contributions. Indeed Kapellasite has about 27% missing 
spins (Cu atoms replaced by Zn ones) on the kagome sites 
and also a few percent of Cu atoms at the center of the 
hexagons. This concentration of Zn on the kagome sites is 
not enough to kill the cuboc2 correlations, the threshold 
being at about 40%^. But triangles with a Zn substitu- 
tion (about 42%) may favor the formation of singlets and 
thus freeze entropy at low temperature. (The spins at the 
center of the hexagons in a symmetric environment ex- 
perience a magnetic field at second order in perturbation 
and are probably only weakly coupled to their neighbors.) 
The quantum description of such a phenomenon is for 
the moment out of the possibilities of any approach and 
theoreticians have to wait for the synthesis of compounds 
with less defects. On the other hand, taking into account 
these two types of defects would probably allow to fit the 
magnetic susceptibility to lower temperatures. (Indeed, 
the pade approximants of the susceptibility converge to 
lower temperature than T m i n ~ 20if , the temperature 
where the pades and the experimental \ data points do 
not agree anymore.) 

A last limit on this determination is the neglect of 
Dzyaloshinskii-Moriya interactions. In fact the lack 
of inversion center on the magnetic bounds allows for 
Dzyaloshinskii-Moriya interaction of spin orbit origin. In 
the cuprates in general these couplings are estimated 
of the order of 1/10 of the super exchange couplings 
and in Herbersmithite they were measured of the or- 
der of 0.06(2)1<PI It is well known that the influ- 
ence of this s mall c oupling is of crucial importance in 
Herbertsmithite^E^, in which it very plausibly explains 
the critical behavior of the compoun d 16 * 17 * 19 ! But in 
that case the influence of this tiny coupling is indeed em- 
phasized by the fact that the anti-ferromagnetic Heisen- 
berg hamiltonian is classically at a multi-critical point 
and quantum mechanically nearby a quantum critical 
point. 16 The situation in Kapellasite is quite different: 
whereas neutron scattering in Herbersmithite is essen- 
tially featureless^^ the experimental evidence of short 
range cuboc2 correlations in Kapellasite is strong^, and 
the results of the present analysis (see Fig 10 ) point in the 



unambiguously but all the fits agree to locate the sys- 
tem in the cuboc2 range, very far away from any critical 
point (the cuboc2-ferro transition is a strong first order 
transition^). Extending the present fit to take into ac- 
count Dzyaloshinskii-Moriya interactions would slightly 
change the exchange parameters but not move the sys- 
tem away from the present phase. With these caveats in 
mind, we think that the present model is the best effec- 
tive model that we are able to build. 



same direction: the J2 and Jd parameters are not fitted 
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Appendix A: HT-series expansion 

We give here the susceptibility HT-series polynomials 
as defined by Eq(3j 



(>2,^3) = 1 

Pi (^2,^3) = -2 - 2z/ 2 - "3 



P2 (>2,^3) 
P3 (>2,^3) 



I Z/3 + 16 Z/ 2 + 8 Z/ 2 Z/ 3 + 4 z/ 2 



(Al) 
(A2) 
(A3) 



-3-48 z/ 3 - 81 v 2 - 3 z/f - 48 z/f z/ 3 + 2 z/f 



12 z/ 2 z/f - 114 v 2 v^ - 96 z/f - 12 z/f (A4) 
P4O2, ^3) = -4 + 176 z/ 3 + 172 1/2 — 4 1/| — 32 z/ 2 z/f 
+216 z/f z/ 3 + 5 z/f + 260 z/f z/f + 396 z/ 2 z/f + 1188 1/31/3 
+432 z/f - 32 z/f + 1094 z/f + 1 108 v 2 z/ 3 + 260 z/f ( A5) 



^5(^2,^3) 



-202 + 225 z/ 3 + 605 v 2 - 680 z/f z/ 3 - 202 z/f 



21 z/f - 270 z/f z/f + 60 v 2 v\ - 3245 z/f z/f - 11430 z/f z/ 3 



~r 60 z/ 3 7970 z/f z/f - 410 z/ 2 z/f - 1360 z/f - 270 z/f 

- 8465 z/ 2 z/f - 17230 z/f z/ 3 - 11645 z/f - 2935 z/f 

- 7630 z/ 2 z/ 3 - 5595 z/f (A6) 
P6O2, ^3) = 1513 + 4104 v 2 - 4206 z/ 3 - 14253 z/f 

+ 122658 z/ 2 z/f + 187998 z/f z/ 3 + 15912 z/ 2 z/ 3 + 132562 z/f 
+ 13170 z/f + 20210 z/f + 98088 z/f + 251238 z/f z/ 3 
+ 40632 z/ 2 z/f + 183249 z/f z/f - 2646 z/ 2 z/f + 104100 z/f z/ 3 
+ 132630 z/f z/f + 38550 z/f z/f - 4170 z/f z/f + 18080 z/f z/f 



26370 z/f z/f + 3384 z/f z/ 3 + 1104 z/ 2 z/f - 4170 z/f 



1104 z4 +6768 z/f + 1513 z 



399 



(A7) 



P7O2, ^3) = 13844 - 151620 v 2 - 74704 z/ 3 + 139083 z/f 



1966153 
2 

189371 



z/ 2 z/f - 1293383 z/f z/ 3 + 299964 z/ 2 z/ 3 



139755 z/? - 371574 v\ 



4759391 



o o 6770421 9 9 
3589292 z/f z/ 3 - 1154895 z/ 2 z/f z^f 



82789 



z/ 2 z/f - 3827614 z/fz/ 3 - 3512439 z/f z/; 



3 7 ,2 

3 



1749132 z/f z/f - 86618 z/f z/f - 990822 z/f z/f 



1991059 z/f z/f - 775061 z/f z/ 3 + 7651 nu 2 v\ 



354109 z/f z/f + 938 v 2 v\ 



67319 



' 2 M 



53018 1/? 



'2^3 



3738 z/f z/f - 262773 z/f z/f + 13615 z/f - 3738 z/f 



1275309 



z/ 9 5 - 106036 z/? + 938 v% + 13844 z/J + 160 vl 



v 2 -r iuu ^ 3 
(A8) 



5(^2,^3) = -186286 - 137536 v 2 + 1145568 1/3 
+ 6324260 z/f - 2088608 z/ 2 z/f - 4756328 z/f z/ 3 
+ 857352 z/ 2 z/ 3 - 10846320 z/f - 1664724 z/f 
+ 1483992 z/f + 16142374 z/f + 41699268 z/f z/ 3 
+ 19842784 z/ 2 z/f + 40719726 z/f z/f 



+ 6650764 z/ 2 z/f + 67176864 z/f z/ 3 + 82220052 z/f z/f 

+ 42915800 z/f z/f + 10425242 z/f z/f + 46137688 z/f z/f 

+ 65417170 z/f z/f + 54159980 z/f z/ 3 - 1232896 z/ 2 z/f 

+ 18554568 z/f z/f + 127572 z/ 2 z/f + 6383884 z/f z/f 

+ 3657208 z/f z/ 3 - 1281328 z/f z/f + 26510312 z/f z/f 

+ 2892272 z/f - 870832 z/f + 32058256 z/f + 4069692 z/f 

+ 248420 z/f + 764352 z/f - 61008 z/f - 186286 z/f 

+ 11421 z/f + 248420 z/f z/f - 592716 z/f z/f 

+ 1625812 z/f z/f + 4060012 z/f z/f + 811044 z/f z/f 

- 61008 v 2 v\ + 382176 z/|z/ 3 (A9) 
P9O2, ^3) = -2329677 + 26960814 v 2 + 11526543 z/ 3 

- 56758545 z/f + 128462472 z/ 2 z/f + 154668708 z/f z/ 3 

- 111054150 z/ 2 z/ 3 - 115342752 z/f - 36200385 z/f 
+ 64998198 z/f + 319743243 z/f - 220814460 z/f z/ 3 

- 147995460 z/ 2 z/f - 152834616 z/f z/f - 244445922 z/ 2 z/f 
822646359 z/f z/ 3 - 1365041772 z/f z/f - 706484844 z/f 

._ o o 



2 /y 3 



- 436794300 z/f z/f - 1401175632 z/f z/f 

- 1793770353 z/f z/f - 1427828508 z/f z/ 3 + 3083328 z/ 2 z/f 

- 944463996 z/f z/f + 2878929 z/ 2 z/f - 456663087 z/f z/f 

- 576786222 z/f z/ 3 + 24409008 z/f z/f - 1226646711 z/f z/f 

- 81429093 z/f + 14618331 z/f - 580389138 z/f 

- 314815344 z/f - 3893763 z/f - 55709883 z/f 

+ 1180962 z/f + 7338168 z/f - 311904 z/f + 7769322 z/f z/f 

- 13919346 z/f z/f - 172316772 z/f z/f - 380586006 z/f z/f 

- 287537220 z/f z/f - 965718 v 2 v\ - 27272862 z/f z/ 3 
+ 1180962 z/f z/f - 5141667 z/f z/f + 9044406 z/f z/f 

- 59040405 z/f z/f - 33586860 z/f z/f - 20688669 z/f z/f 

3669084 z/f z/ 3 - 2329677 z/f + 37370 z/f 

(A10) 



-311904 z/ 2 z/f 



)z/f 



P10O2, 0) = 44494564 + 32699900 v 2 - 1607336300 . z 
+ 4682885400 z/f - 1969984450 z/f - 3722864284 z/f 
+ 12819641560 z/f + 2641862210 z/f + 853326455 z/f 

- 109501560 z/f + 44494564 z/f (All) 
Pio(0, ^3) = 44494564 - 323940580 z/ 3 + 753910650 z/f 

- 746061580 z/f + 248953155 z/f + 673343648 z/f 

- 251041900 z/f + 86401210 z/f - 21228160 z/f 

+ 4486600 z/f - v f (A12) 

\ 5959595279 19215984161 2 
Pn(0, 1/3) = 568071766 z/ 3 + z/f 

49178193933 o 4 , 
z/f + 33878328495 z/f - 30148742943 z/f 

fi 7 1749626615 « 
+ 8479066530 z/f - 2914976526 z/f + z/f 

- 246955709 z/f + 58891206 z/f° - 8569254 z/f 1 

(A13) 
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We give here the specific heat HT-series polynomials as 
defined by Eqjl2} 



90(^3,^3) = 
91(^3,^3) = 

92(^3,^3) =3 + 3i/| + ^1/3 



9 3 

■4 



(A14) 
(A15) 

(A16) 
(A17) 



9 

<?3(>3,^3) = -27^2 -54^3^+2 
44 (>3, ^3) = -153 + 144 i/ 2 + 108 i/ 3 + 144 v 3 v 2 - 90 v\ 

o /i 00 45 

^3 z/ 2 — 153 z/J — 90 ^3 ^2 — ^~ ^ : 



- 252 1/3 i/f - 72 v\ v 2 



(A18) 

45O3, ^3) = 3300 v 2 - 900 v 3 - 2850 i/f + 5700 v 3 "2 

- 150 v\ + 1650 z/f - 4650 1/3 i/| + 750 v\ v 2 - 450 z/f 
+ 6000 v 3 v\ - 1350 v\ v\ + 3300 z/f i/ 2 - 300 z/f z/f 

- 450 z/f z/| - 225 z/f (A19) 
^eOs, ^3) = 32085/2 - 24570 z/ 2 - 20655 i/ 3 - 21735 z/| 

+ 1620 ^3 ^2 + + 19890 z/f - 131085 1/3 z/f 

+ 43200 z/f ^ 2 - 9225 z/f + 9450 v\ + 46980 ^3 z/f 

- 81675 z/f z/f - 2430 z/f v 2 + v\ - 32805 i/ 3 v\ 

+ 54000 z/f z/f - 28080 z/f z/f + 19710 z/f v 2 + — ^— z/f 

+ 16875 z, 2 v\ + 4950 z/f z/f + ^ i/f z, 2 + ^ z^ 6 

(A20) 

$7(^3, ^3) = -10143 - 644301 v 2 + 269892 i/ 3 + 1052226 z/f 

- 934920 ^ 3 ^ 2 - 75411 z/f - 460404 z/f + 1893654 i/ 3 z/f 

- 31752 z/f ^ 2 + 122598 z/f + 119511 z/f - 1000629 z/ 3 

o 9 o 4 930069 * 
+ 242991 z/f v\ - 1032381 z/f v 2 - 31311 z/f — z/f 

+ 269010 i/ 3 f 2 + 827757 v\ v\ - 257544 z/f z/f 

— v\ v 2 + 58653 z/f - 1421343 i/ 3 v\ 

52479 

- 187425 z/f v\ - 980343 z/f z/f + — ^— v\ z/f 

- 227997 z/f v 2 - 10143 z/f + 7938 z/f z/f - 38808 z/f z/f 

57771 7 



-66150z/fz/f + 58653 z/f z/f , 



^ ^3 ^2 t ouuu^» ^3 ^ 2 -r ^ ^3 

s(^3, ^3) = -2859213 + 6178704 z/ 2 + 5449416 v 3 
+ 9429168 z/| - 9235296 1/3 i/ 2 - 4847052 v\ 
- 20274912 z/| + 53705904 z/ 3 i/| - 21687792 v\ v 2 

- 5012952 vl + 9225132 z/| - 62253744 u 3 v\ 

- 48280344 v\ v\ + 3022992 z/| v 2 - 1229970 v\ 



(A21) 



1 ^^^^ -3 -^^^^.^^3 

- 1534176 z/| + 56331912 v 3 v\ - 65766624 v\ v\ 
- 40230624 v\ v\ - 12512304 v\ v 2 + 1025136 v\ 



- 125076 , 



- 626304 z/| - 7418544 v 3 v\ + 43961232 z/ 2 i/| 

- 8318688 z/| i/| + 23324616 v\ v\ - 130368 v\ v 2 

- 125076 v% + 10753344 z/ 3 i/| - 12882912 z/ 2 i/| 
+ 10512768 v\ v\ - 13091568 v\ v\ + 4212600 1/| v\ 

- 3452400 v% v 2 - 2859213 v\ - 3844932 v\ v% 

- 741888 v\ v\ - 2292990 v\ v\ - 654192 v\ v\ 

• 'T' -f (A22) 
^9(^3, ^3) = 5600664 + 184435056 1/ 2 - 100067400 z/ 3 

- 433650996 v\ + 227967048 1/3 i/ 2 + 78952644 z/ 2 

+ 179433792 v\ - 607203540 v 3 v\ - 42816600 v\ v 2 

- 63694188 v\ + 95703120 v\ - 144709416 v 3 v\ 

+ 485374032 v\ v\ + 276310440 1/| i/ 2 + 35501328 v\ 
+ 123924168 z/| + 441623340 z/ 3 ^2 - 1670214816 v\ v\ 
+ 372083544 1/| z/| - 127644984 1/| z/ 2 - 24001596 v\ 

- 70144056 v% + 144516960 z/ 3 i/| + 1235200104 v\ v\ 

- 631218096 v\ v\ + 54034452 v\ v\ + 115234488 1/| z/ 2 
+ 10998828 z/| + 146195604 v\ - 127553616 v 3 v% 

- 684430884 v\ v\ + 852307920 v\ v\ - 593753004 v\ v\ 
+ 250454268 v\ v\ + 4221396 1/| i/ 2 - 8027748 1/| 

+ 442063656 v 3 v\ + 75600216 z/| i/| + 423090216 1/| v\ 
+ 106291116 z/| z/| + 101241576 1/| i/| + 22410432 v% v\ 
+ 12159720 z/J i/ 2 + 5600664 z/| + 4667544 v\ v\ 
+ 32428512 z/| i/| - 21895272 v\ v\ + 2416392 1/| i/| 

- 17600328 1/| z/| - 8027748 v\ v\ - 1161810 v\ 

(A23) 

(7io(z^2,0) = 559Q ^ 5695 _ 2276071650 v 2 - 4567367250 v\ 

+ 14999658300 1/| - 11395200825 v\ + 3512454300 z/| 

- 4851077175 1/| + 1271847150 v\ - 191251800 v\ 

(A24) 



1559095695 ln 

7/„ 



^ >. 559095695 _^ ooooor7r 3707972775 2 
9io(0, ^3) = g 1968883875 z/ 3 + = vl 



- 2766779100 z/| + 953949150 1/| - 929445975 1/| 

a 7 59072625 
+ 73620225 z/| - 111223125 i/J — 



45148455 10 
: ^3 



4 

9ii (0, z^ 3 ) = 3188690010 + 49046431005 v 3 

127470667995 2 103764854985 3 
7> 4 + « ^1 



(A25) 



2 

I - 15771729270 v% 



z z 

- 41081190975 z/| + 16829360625 1/| _ . . . . _ , . 

10254468675 7 « 
z/J - 3039511530 z/| + 1205526630 u 



390259485 n 

: 



(A26) 
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Appendix B: Physical Pade Approximant (PPA) 

For fixed values of the coupling constants Ji, J2 and 
Jd, we evaluate the order- n HT-polynomial of the mag- 
netic susceptibility, A^ T (/3), or of the specific heat, 
Cy™' (e), where e is the energy per spin. From a poly- 
nomial f n (x) of order n we calculate n + 1 rational frac- 
tions p n -d(aO/(Zd(#)] with d running from to n. A PPA 
must have no zero in p and q in the whole interval of 
variation of x, i.e. in [0, 00] for X^ T (/3) and [eo,0] for 
Cy^ n ' HT (e), where eo is the ground state energy. Solu- 
tions with zero in the range of interest are discarded. By 
varying the coupling parameters, the number of PPAs 
may eventually change. Thus all functions build on the 
sum over the PPAs are discontinuous. Unfortunately this 
prevents from using minimization powerful method. This 
is a price to pay when using PPAs. 



Appendix C: Determination of the parameters A 
and B of Qx 



where X means the average value over the set of temper- 
atures. Then, the best PPA (highest Qx,ppa) leads to 
the first values (A PPA , B PPA ) which are used to compute 
the Q x ,ppa and Eqs |C3|[C4| are solved as: 

(C8) 
(C9) 

A = (TF^ 2 ) - (t*) (^p7) (CIO) 

where (X) means the average value over the set 
of temperatures and over the PPA with the weights 

Q^ppa/(1 — Q^ppa) 1-2 ^- The procedure is iterated 
until convergence by calculating the new weights at the 
new A and B. The convergence is quick and a couple of 
iterations are sufficient for a relative precision of 10 -5 on 
A and B. 



From Eq|6]or Eq|l4j we have: 

QX = 5. l + \Rx.PPA/e\ r 



PPA 



(CI) 



Rx,ppa = [AF PPA (T k ) + BT k - F fc exp ] 2 (C2) 



where X = X or Cy and B is for Cy. 



dQ 



x 



dQ 



X,PPA 



dA ^ OA 

{PPA} 



2e2NT {lk } (I-Qx.ppa) 1 " 2 /'' 



.1+2/r 
^X,PPA 



J2 l(AF PPA (T k ) - BT k - F fc exp ) F PPA (T k )] 

(C3) 



k=i 



an n 1+2 / r 

dB _ 2^7V T{ A 1} (l-Q X5 p PA )i-^ 



53 [(AF PPA (T k ) - BT k ~ F fe exP ) Tk\ (C4) 



fc=l 



where F stands for XTjC or CV/T. We look for A and 
5 that cancel out these derivatives. First, A and £> are 
evaluated for each PPA: 

Appa =^ (tF-pTF P pa - T^Fppa F«p) (C5) 
^ppa =^ (FppaTFppaF-p - TF^F^j (C6) 
A=TTVp7 2 -T^^|p7 (C7) 



Appendix D: Pade approximant for Cy 

We first determine the HT-series of s(e) the entropy 
per spin versus the energy per spin 6 : From the HT-series 
expansion of Cy(T) = ^Z^ =2 a iP l ^ with j3 = 1/T, we 
obtain the HT-series of s(T) and e(T) as: 

00 spin /rpi 



, c s r(V) 



poo 

s(T)=ln2-y dT' 



= ln2 - ^ -ft + (Dl) 



2=2 

e(T)=- J dT'cf n (T') 

n-l 



i=2 



(D2) 



where we use s(T = 00) = In 2 and e(T = 00) = 0. The 
HT-series expansion of s(e): 



s(e) = ^bie* 



(D3) 



i=0 



is obtained order by order using (Mapple style): 



for i from to n do 
ss :=0; 

for j from to i-1 do 

ss :=ss+b [j] *coef f (E~j ,beta, i) ; 
od; 

b[i] :=(coeff (S, beta, i) -ss) /coef f (E~i , beta, i) ; 
od; 

where S and E and stand for the HT-series of s(T) and 
e(T). 
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The low-T leading term of Cy{T) is assumed to be 
(T/To) a . Then s(e) oc (e — eo) 1 ^ for e around eo where 
eo is the ground state energy and \i = 1 + 1/a. We define 
a function: 



G(e) 



s{eY 



e - eo 

The HT-series expansion for G(e) is obtained from: 
(In 2)^ 



G(e) 



P(e) = 



e 



7.1 



=0 



i=0 



In 2 



i=2 



(D4) 

(D5) 
(D6) 



where 6* = ^/ln2 (b t = 0) and = r(/i + l)/r(/i + 

1 — i) = — l)...(/i + 1 — i). Keeping only terms up to 
order n defines G HT (e). Note that P(e) 2 starts at order 
2i. 

Then G HT (e) is transformed in all possible pade ap- 
proximants noted G^ T (e) where n — d and d are the 
numerator and denominator degrees. We keep only 



the physical pade approximants (PPA) denoted G^* T (e) 
whose numerator or denominator have no zero inside the 
[eo,0]. The value G^* T (eo) is related to cq (see EqpTj): 



G£ T (eo) = c; 



0,d* 



-iCt + 1 



(D7) 



From a G^* T = N n _d* ( e )/Dd* (e), we obtain first s(e) 
and its first derivatives and then the specific heat: 



Sd*(e) 
s d *(e) 
Sd*(e) 



(e - eO) 
1 



7V(e) 



-eO 

<§d*(e) 



Die). 
N'(e) _ D'(e) 
N(e) D(e) 

N"{e) D"{e) 
+ N(e) ~ D(e) 



\N'(ey 


2 


ri)'(e)- 


1 




+ 




(e-e ) 2 



(D8) 
(D9) 

(D10) 



Then we deduce /3(e) 
-K,(e)] 2 /s' d ',(e). 



1/T(e) = 4(e) and Cy in (e) 
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